clear all
%clc

nMC = [1:250];

theta_all = [];
ratio_trimmed_mean_all = [];
ratio_trimmed_mean_all_first = [];
mean_ratio_all = [];


for i=nMC
    
    tempp = num2str(i);
    temp = strcat('MC_',tempp);
    temp2 = strcat(temp,'.mat');
    
    load(temp2)
    
    theta_all = [theta_all;
        theta_NPD'];
    
    ratio_trimmed_mean_all = [ratio_trimmed_mean_all;
                              ratio_trimmed_mean_temp];

    ratio_trimmed_mean_all_first = [ratio_trimmed_mean_all_first;
                              ratio_trimmed_mean_temp_first];

end

theta_NPD = theta_all';
%save('results.mat')

display('*** RATIO OF TRIMMED MEANS -- SECOND ***')
mean(ratio_trimmed_mean_all)
std(ratio_trimmed_mean_all)

display('*** RATIO OF TRIMMED MEANS -- FIRST ***')
mean(ratio_trimmed_mean_all_first)
std(ratio_trimmed_mean_all_first)